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We propose blue-detuned optical traps that are suitable for trapping of both ground state and 
Rydberg excited atoms. Addition of a background compensation field or suitable choice of the trap 
geometry provides a magic trapping condition for ground and Rydberg atoms at the trap center. 
Deviations from the magic condition at finite temperature are calculated. Designs that achieve less 
than 200 kHz differential trap shift between Cs ground and 125s Rydberg states for 10 ^iK Cs atoms 
are presented. Consideration of the trapping potential and photoionization rates suggest that these 
traps will be useful for quantum information experiments with atomic qubits. 
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I. INTRODUCTION 



While ground state neutral atoms interact 
only weakly due to small van der Waals and 
magnetostatic dipolar interactions several re- 
cent experiments have shown that Rydberg ex- 
citation can be used to turn on strong interac- 
tions suitable for quantum gates and entangle- 
ment generation PHI . Following these develop- 
ments Rydberg mediated quantum gates Q are 
currently being studied intensively as a route to 
scalable quantum information processing [f|. 

Recent quantum gate experiments have used 
Rb atoms that are laser cooled and then trans- 
ferred into red detuned far-off resonance opti- 
cal traps (FORTs). Red detuned traps are ad- 
equate for ground state atoms but they present 
several difficulties for experiments that rely on 
Rydberg excitation. The trapping light pho- 
toionizes Rydberg atoms with typical photoion- 
ization rates for few mK deep traps that can be 
faster than radiative decay rates d-Q, and so 
photoionization presents a limit to the usable 
Rydberg lifetime. Furthermore the differential 
light shift of the Rydberg and ground states re- 
sults in a position dependent Rydberg excita- 
tion energy, unless the atoms are cooled to the 
motional ground state of the trapping potential. 
Variations in the excitation energy impact the 
detuning of pulses used for gate operations, and 



thus degrade gate fidelity. To get around these 
limits the trap light is turned off during the gate 
sequence, and then turned on again afterwards. 
Provided the atoms are sufficiently cold, and 
the Rydberg gate lasts only a few fxs, turning 
the trap on and off docs not lead to appreciable 
heating or atom loss out of the trap. 

In a multi-qubit experiment we envision an 
array of optical traps, each holding a neutral 
atom. In most implementations using lattices, 
or trap arrays generated with diffractive optics, 
it is not possible to control the trap intensity on 
a site by site basis. It would therefore be nec- 
essary to turn off all traps whenever any atom 
is Rydberg excited. This is problematic for im- 
plementations with many qubits and it is there- 
fore of interest to find traps that work for both 
ground and Rydberg state atoms. Since the Ry- 
dberg polarizability is that of a free electron 
and is negative, a stable trap must be a dark 
region surrounded by light, and the trap wave- 
length should be chosen so also the ground state 
polarizability is negative. In the alkali atoms 
this implies tuning to the blue of one or both 
of the first resonance lines @. Although a blue 
detuned trap can be attractive for both ground 
and Rydberg atoms, trap depth matching is still 
an issue due to state dependent differences in 
the magnitude of the polarizability and due to 
the different spatial extent of the Rydberg wave- 
function compared to the ground state atom. 
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For high fidelity quantum gates we expect to 
access Rydberg levels with principal quantum 
number n > 100@ and it is therefore necessary 
to consider the local trapping potential aver- 
aged over the Rydberg electron wavefunctionQ 
which may extend to more than 1 /jm away from 
the nucleus. 

Several authors have considered low fre- 
quency electromagnetic trap designs for ground 
and Rydberg atoms floj| . In this paper we show 
that optical frequency traps can be used for 
both ground and Rydberg state atoms, and that 
position dependent differential light shifts can 
be minimized in what we refer to as "quasi- 
magic" trap geometries. In Sec. [TT] we present 
three alternative designs for blue detuned opti- 
cal traps. In Sec. Mil we calculate the Rydberg 
trapping potential and identify magic trapping 
conditions. Representative numbers are given 
for Cs atoms. Photoionization rates are calcu- 
lated in Sec. IVland we conclude in Sec. EU 



II. BOTTLE BEAM OPTICAL TRAPS 

Wavelength regions where the ground and 
Rydberg state polarizabilities are the same sign 
are to the blue of the first resonance lines in al- 
kali atoms. Calculated polarizability curves for 
the heavy alkalis Rb and Cs are shown in Fig. 
[U The curves for the 50d Rydberg state are 
within a few percent of the value found from 
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the free electron polarizability a e = - e a , ex- 
cept near the resonance with the first excited 
p level. We see that for both elements there is 
a matching wavelength to the red of the sec- 
ond resonance lines at approximately 430 nm 
for Rb and 470 nm for Cs. The ground state 
vector polarizabilities are also very small at this 
wavelength which implies small rates for hyper- 
fine or Zeeman state changing Raman transi- 
tions. Although these wavelengths might there- 
fore appear attractive for trapping ground and 
Rydberg states they are not useful due to the 
need to account for the different spatial extent 
of the ground and Rydberg wavefunctions. As 
we will see in Sec. IIIII it is preferable to work 
at longer wavelengths for which the ground and 



5s-6p 5p 3/2 -50d 5/2 5s-5p 




wavelength [\im) 

FIG. 1. (color online) Scalar polarizability of 
ground and Rydberg states of Rb and Cs (solid 
lines) and vector polarizability of the ground state 
(dashed lines). 



Rydberg polarizabilities are both negative, but 
the ground state polarizability is much larger in 
magnitude than that of the Rydberg state. 

Several methods have been used to produce 
bottle beam traps (BoBs) that have an inten- 
sity null surrounded by light in all directions 
[l2|, We have investigated in detail 

the three configurations shown in Fig. [2] The 
Gaussian interference BoB in Fig. [2^) makes 
use of the interference of two TEMoo Gaussian 
beams with different waist sizes Wi, W2 [f3j . 
The crossed vortex BoB [13] in Fig. [2Jd is formed 
by two Laguerre-Gaussian Lq beams with or- 
thogonal polarizations that cross with an angle 
of 29. We have recently demosntrated trapping 
of single ground state Cs atoms in both of these 
BoB traps [3. 

The third dipole trap [2J; is created by 
four parallel tightly focused TEMoo Gaussian 
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beams. The four beams with a waist size wq are 
spaced on a square with sides d. Each beam has 
the same polarization as its diagonal neighbor 
and has orthogonal polarization to that of its 
nearest neighbors. This polarization configura- 
tion minimizes the effects of interference. Both 
the waist size wq and beam spacing d are on 
the \im scale. The overlap of the four beams 
forms a potential barrier around the center of 
the square in the x — y plane. Diffractive spread- 
ing of the Gaussians also creates a trapping po- 
tential along z, thus forming a 3D BoB trap. 
This latter configuration is of particular inter- 
est for forming tightly packed BoB arrays. 

For each trapping geometry atom localization 
near the trap center can be quantified by an 
expansion of the potential U(r) — —^-^al(r) 
near the trap center. Here a is the scalar polar- 
izability and / is the intensity at position r. The 
intensity distributions for the different trap con- 
figurations are calculated in Appendix A. For 
the Gaussian interference BoB we find near the 
origin 
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where a is the scalar polarizability of the atom, 
A is the trapping wavelength, and Pi , P2 = 
(W2/W1) 2 Pi are the powers of the beams with 
waists w\ , W2 rcpcctivcly. The trapping poten- 
tial is axially symmetric and quartic in the xy 
plane and quadratic along z. The total trap 
power used in Fig. [2] is P = P\ + Pi- 

For the crossed vortex BoB an expansion 
about the origin yields 
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where P is the total power of the two beams, 
and w is the focused waist size. This trap is 



Telescope 




b) 



spiral 
phase plate 

Calcite 




| Feedback 
Control 



3 



FIG. 2. (color online) Setup of blue detuned dipole 
traps: a) Gaussian interference BqB[I3||. b) crossed 
vortex BoB [l2T] . and c) Gaussian lattice. 



quadratic in all directions. 

For the Gaussian array BoB the expansion 
along x and z is 

2w 2 - d 2 



U(x,0,0) = -U e~^ 1 



x A + 0(x 4 ) 



U{0,0,z) = -u e: 
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with Un = 8aP 2 . Here P is the power of each 
beam in Fig. In an array implementation 
each beam is shared between four neighboring 
trapping sites, so a total power of only P per 
site is requircd(this neglects a small correction 
due to the rows at the edge of the array). 

Trapping potentials of the three dipole trap 
configurations for the ground state of Cs are 
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FIG. 3. (color online) Trapping Depth of blue de- 
tuned dipole traps in the x — y (left column) and 
x — z (right column) planes for Cs 6s, A = 780 nm. 
a = -235. x 10" 24 cm 3 , power P = 50 mW for 
a) Gaussian interference BoB wi = 2 /im, W2 = 
3.78 /xm, b) crossed vortex BoB w = 3 /xm, 6 = 8.6°, 
and c) Gaussian lattice trap w = 1.5 /xm, d = 4 /xm. 



III. PONDEROMOTIVE POTENTIAL 
OF TRAPPED RYDBERG ATOMS 

High n Rydbcrg atoms with n > 100 have 
electron wavefunctions that are comparable in 
spatial extent to the trap potentials shown in 
Fig. O The AC Stark shift of Rydberg atoms 
can therefore no longer be approximated by 
U = — 2^~7 a ^i with / the local intensity at the 
nucleus. We need to consider the ponderomo- 
tive energy of Rydberg atoms in a field of vary- 
ing intensity. The pondcromotivc shift is the 
time averaged kinetic energy of a free electron 
in an oscillating electric field. For a field of the 
form Ecos(ujt), the ponderomotive energy is 



e 2 \E\ 2 
Am e u) 2 



where — e and m e are the electron charge and 
mass respectively. Using / = ^|i?| 2 where c is 
the speed of light we can write the ponderomo- 
tive energy of a free electron as 
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Then the Hamiltonian of a Rydberg atom in an 
oscillating electromagnetic field can be written 
as 

H F + U P (R + r)V(r; R) = E R {R)^{r; R), 



plotted in Fig. [3] and the trap oscillation fre- 
quencies along different axes are listed in Table 
HI We see that all three designs provide trans- 
verse oscillation frequencies of a few tens of kHz 
and longitudinal oscillation frequencies of a few 
kHz. The vortex and Gaussian lattice traps re- 
sult in quite similar frequencies and trapping 
depth. The Gaussian interference BoB is about 
3 x shallower for the same optical power and has 
the poorest axial confinement. 



where R is the center of mass coordinate of the 
atom, and r is the coordinate of the electron 
relative to the center of mass. Using first order 
perturbation theory, and supposing there is no 
degeneracy involved, the energy shift of a Ryd- 
berg atom in state j is [t| 

A^-(R) = J d 3 rUp(-R + r)\^(v;R)\ 2 

= 9 ^ 2 / d 3 rI(R + r)|^(r; R]Qi) 

This expression is valid provided the pondero- 
motive potential varies over distance scales that 
are larger than the wavelength of the Rydberg 
electron. This is well satisfied for the potentials 
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TABLE I. Oscillation frequencies and the trap potential at the lowest saddle point for the BoB traps, each 
with the same total power of 50 mW, with trap parameters from Fig[3] 



design 


w x /27T (kHz) 


ujy/2-n (kHz) 


u) z /2n (kHz) 


U/k B (A*K) 


Gaussian 
Interfer- 
ence 


62.5^ 


62.5^ 


0.315 


60 


crossed 
Vortex 


29.4 


29.8 


4.42 


225 


Gaussian 
Lattice 


15.4 


15.4 


2.79 


256 



a The Gaussian interference BoB is anharmonic in the radial direction. The vibration frequency was calculated by 
setting the particle energy to 1/10 of the trapping potential. 



we consider. At n = 150 the electron wave- 
length is about 50 nm which is less than 10 % 
of the wavelength of the light creating the trap. 
In addition it is necessary that the pondcromo- 
tive shift is everywhere small compared to the 
energy spacing of Rydbcrg levels. For the 150s 
state, which is the highest we consider in the ex- 
amples below, the closest state is 1 46/7/2 which 
is 1.6 mK away. Looking at Fig. 2] the largest 
ponderomotive energy seen by a 150s atom for 
the traps we are considering is about 200 /iK. 
The ratio of energy scales would thus imply a 
higher order correction ~ 0.12. 

In fact this naive estimate is overly pes- 
simistic since the coupling between ns and (n — 
4)/ is strongly suppressed by the trap geome- 
try. The dipole traps in Fig. [5] all have spa- 
tial reflection symmetry so that the coupling 
matrix elements between the ns and (n — 4)/ 
Rydberg states are exactly zero when the atom 
is at the origin. For the low atom tempera- 
tures expected for Cs, the atom will be very 
near the center of the trap and thus this cou- 
pling will be strongly suppressed. The next clos- 
est states are in the (n-4) degenerate manifold 
(I = 4, 5, ...) which are separated from the 150s 
state by about 5 mK which is more than a fac- 
tor of 25 larger than the light shift. This will 
lead to a second order perturbative correction to 
the energy shift of approximately 1 part in 25. 
Since the Rydberg level spacing scales as 1/n 3 , 
with n the principal quantum number, the error 
will be even smaller for lower levels. 



We will calculate the wavefunctions x/j® using 
a model pseudo potential method. The poten- 
tial form adopted here is [l5[ 

Zi(r) 



Vi(r) 



2r 4 



-(r/r a y 
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where Z t (r) = 1 + 36e~ Q i r 
a c i = 15.81, r c = 2.0, and all the other pa- 
rameters are listed in Table [XXJ To verify our 



TABLE II. Parameters for the Cs model potential 



I 





1 


2 


3 


4+ 




3.49625 


3.73801 


3.45092 


3.43592 


3.43592 




9.57499 


9.56664 


9.52285 


9.54285 


9.54285 




1.41409 


1.34016 


1.58147 


1.62147 


1.62147 



calculation of the wavefunctions, we reproduced 
the planewave photoionization cross sections 
listed in Q. The fine structure of Cs adds less 
than 0.1% correction to the ponderomotive 
energy shifts, so we ignore fine structure 
corrections in this paper. 

Figure |4] gives sample calculation results for 
ns Rydberg levels with n = 100, 125, 150. We 
see that as n increases the effective trapping 
potential gets smaller and smaller. This is be- 
cause the large electron wavefunction averages 
over the intensity distribution of the trap ac- 
cording to Eq. ((4]) which washes out the poten- 
tial minimum. If the trap parameters are not 
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FIG. 4. (color online) Potential energy of Cs ground 
and ns Rydberg states with n = 100, 125, 150 in a) 
Gaussian interference BoB, b) vortex BoB, and c) 
Gaussian lattice BoB. Trap parameters the same as 
in Fig. [31 



chosen correctly, as is the case in Fig. @J>), the 
trap could be repulsive for high n even though 
a e is negative. Even when the trap is attrac- 
tive for Rydberg states the ground to Rydberg 
trap shift for an atom at R = is not negligi- 
ble. This shift increases with n and is propor- 
tional to the light intensity. In an experiment 
with an array of traps this would imply that 
the Rydberg excitation energy would vary from 
site to site due to intensity variations across 
the array. To minimize this effect we seek trap 
parameters for which the R = trap induced 
shift vanishes. We will refer to this in what fol- 
lows as "quasi-magic" trapping. A quasi-magic 
trap will give an intensity independent excita- 
tion shift for atoms at the trap center (or for 
atoms in the motional ground state with slightly 
different compensation parameters) and only a 
small shift for sufficiently cold atoms. We quan- 



tify the notion of small in the following section. 

IV. MAGIC CONDITION FOR ZERO 
TEMPERATURE ATOMS 

Inspection of Fig. [1] shows that apart from 
wavelengths that are very close to the second 
resonance lines the magnitude of the ground 
state polarizability is larger than that of the 
Rydberg state. Conversely Fig. @] shows that 
the trapping potential at R = is larger for 
Rydberg states than for ground states. This 
implies that we can balance the R = trap 
shifts by adding a constant background inten- 
sity that will shift the ground state potentials 
more than the Rydberg state potentials. With 
the correct background intensity / m the differ- 
ential shift will vanish. This is the quasi-magic 
trapping condition. Note that if we were to use 
the wavelengths in Fig. [T] where the ground 
and Rydberg polarizabilitics arc equal we would 
have to add a relatively large background inten- 
sity. At A = 780 nm the ground state polariz- 
ability a is about 5.4 x larger than that of the 
Rydberg state a e which reduces the power re- 
quirement for the background beam by this fac- 
tor. It is possible to work even closer to the first 
resonance line where a/a e is even larger, but 
decoherence rates associated with photon scat- 
tering and differential hyperfme shifts H, in- 
crease correspondingly. We have therefore cho- 
sen 780 nm for Cs as a viable working wave- 
length. 

Using the ground state light shift 
AU ff = --^[/BoB(R) + / m (R)], 

Ze c 

and the Rydberg state shift 



2eocm e ui 2 

J d 3 r[I BoB (R + r) + J m (R + r)]|^°(r; R)| 2 

the quasi-magic condition is simply AU g = 
AUr. Figure [5] shows an example of such a 
magic condition for the crossed vortex BoB. 
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FIG. 5. (color online) Energy shift compensation 
for the crossed vortex BoB with a planewave of in- 
tensity I m = 128 nW/nm 2 , w = 3 fim, 9 = 8.6°, 
and P = 50 mW. 
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Although the additional power required for 
matching is small for a single site, the addi- 
tional light requirement becomes substantial if 
we consider a 100 x 100 or 1000 x 1000 ^m 2 
array. The Gaussian lattice design presents an 
interesting alternative since the light intensity 
is naturally not zero at the trap center. The 
R = intensity changes as we vary the waist 
size or separation of the beams, and by judi- 
cious choice of parameters we can achieve the 
matching condition without adding any addi- 
tional plane wave. Note that the compensating 
intensity is in this case not uniform but is spa- 
tially varying. Figure [5] shows such a self magic 
condition for n = 125. 

For a ground state atom with a low tempera- 
ture, we can estimate the average trap induced 
shift between ground and Rydbcrg states by 
< dU >= |E l=w ^(0,0,0)(rf), where 
the mean square position of the atom found 
from the Virial theorem is < r 2 >= 2 g B Jj , 

dUu = dii(Un —Ug) and U g , Ur are the ground 
and Rydberg state trapping potentials. Fig- 
ure [7] shows that the transition shift decreases 
nearly linearly with decreasing atom tempera- 
ture. This shift would be below 0.2 MHz for 
an atom temperature of 10 [iK which is readily 
achieved using polarization gradient cooling of 
Cs. 



FIG. 6. (color online) Trapping potential (top row) 
and shift difference (bottom row) between Cs 6s and 
125s for a self-magic Gaussian lattice trap with A = 
780 nm, d — 4 /^m, w = 1.57 fj,m, and P — 50 mW. 
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FIG. 7. (color online) Average transition shift be- 
tween ground and 125s states of Cs in Gaussian lat- 
tice trap with d = 4 /im, w = 1.57 /im, P = 50 mW, 
and E/trap = fcs x 300 /jK. 



V. PHOTOIONIZATION RATE 

In this section we calculate the photoioniza- 
tion rate of Rydberg atoms in a BoB trap. Since 
the Rydbcrg electron is not exposed to a uni- 
form intensity field the photoionization calcula- 
tion must be modified accordingly. The local 
photoionization rate R scales as R = cr^j with 
u> the photon frequency. The cross section is 
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a = 2tt' 



dfif 



m1c 2 a dE 



a is the fine structure constant and the deriva- 
tive of the oscillator strength with respect to 
the energy is 



27). 



dfil _ 
dE muj 21, 



1 



TOi l f ,m f 



K/|a-pN)| 2 . (6) 



where the initial state |i) = \rii,li,mi) is a Ry- 
dberg state with principle quantum number rii, 
and the final state |/} = \Ef,lf,rrif) is a con- 
tinuum state with energy Ef = E^ + En, with 
E u = Tnu the photon energy. The magnitude 
of A is normalized to unit peak intensity. Even 
though kr ~ kn 2 ao 3> 1 is large for our pa- 
rameters we may ignore high powers of r in the 
expansion of A when calculating the matrix el- 
ement since, even though the electron's wave- 
function is comparable in size to the photon's 
wavelength, the pho ton absorption takes place 
near the nucleus jl8|. The quadrupole term is 
included due to the small electric dipole transi- 
tion rate for s state atoms near the center of a 
dark trap. We tested the matrix element calcu- 
lation using the full multipolc operator and ver- 
ified that only the dipole and quadrupole terms 
gave a substantial contribution. 

For a planewave field polarized in the x di- 



rection A = 
imatcd by 



e x, and Eq. (J6J) can be approx- 



dftf 2h 1 
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Eq. ([7]) becomes 
dfi f 2mu> 1 



dE 



h 21 + 1 



a ih x 2 , i*k>n , ik z | ,, 
\x + —x + — xy+ — xz\i) 



-EE 

m» l f , mf 

For a spatially inhomogeneous field like 
the Gaussian lattice BoB we decompose into 
plancwaves as 

1 



A(r 



(2tt) 3 
,9k = J d 3 rA(r)e 



d 3 kg k e* kr , 



— ik-: 



The oscillator strength derivative can then be 
written as 

dfi f 2mw 



dE (27r) 6 n 21 



m, l f ,m f 



d 3 kg k [{x) + ^(xx) 



(xy) 



ikz 
2 



{xz} 



To evaluate the matrix elements the wavefunc- 
tions are calculated with the same method as in 
Section IIII1 The radial part of the continuum 
state is normalized to 



' 2m 



sin(fc e r + S), as r 



dE moj 2li + 1 

x 2J zJ \(f\P* + ikxXPx + ik y yp x + ik z zp x \ 

rrii lf,rtif 

Using the following relations 

xp x = (Hxx — xxH) , 

yp* = W( Hx v - x v H ) - h l * 

zp x = ^{Hxz — xzH) + ^l y . 

and dropping the magnetic dipole terms, which 
give no contribution to the photoionization rate, 



where k e is the wavenumber of the free electron 
and 5 is the continuum state phase shift. 

Figure [5] shows the photoionization rate for a 
125s state Cs atom in a Gaussian lattice trap 
which satisfies the self magic condition of Fig. [6j 
The quadrupole term gives less than 3% correc- 
jtipn to the final result shown in the plot due to 
the non-zero light intensity at the trap center. 
The photoionization rate is also substantially 
smaller than the room temperature radiative 
decay rate of the same Rydberg state which is 
about 1800 s _1 . Nevertheless the room temper- 
ature photoionization rate of the 125s state is 
dominated by blackbody radiation which gives 
a rate of about[l]| 20 s _1 at 300 K. Only at 
cryogenic temperatures less than 10 K does the 
trap light induced rate at trap center seen in 
Fig. [5] dominate over the blackbody rate. 
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FIG. 8. (color online) Photoionization rate for 125s 
Cs in a 780 nm self-magic Gaussian lattice dipole 
trap, w = 1.57 /jm, d = 4 //,m, P — 50 mW, and 
f/trap = k B x 300 /iK. 

VI. DISCUSSION AND CONCLUSIONS 

In summary we have presented three designs 
for blue detuned dipole traps that arc capa- 
ble of trapping both ground and Rydbcrg state 
atoms. Using visible or near-infrared trapping 
wavelengths, and alkali atoms with tempera- 
tures < 100 /J.K, these traps are capable of /mi 
scale atomic localization in three dimensions. 
We have calculated the pondcromotive poten- 
tial energy of trapped Rydberg atoms, the im- 
portance of which has been demonstrated in re- 
cent experiments [2(| , and shown that it is pos- 
sible to match the ground and Rydberg state 
trap shifts for atoms at the center of the trap. 

One attractive feature of these optical traps 
is that they can be replicated easily in two di- 
mensions with a diffractive beam splitter. In 
this way the traps could be used in experiments 
that require control over individual sites of a 
closely spaced two-dimensional atomic array as 
in[2ll|. This type of holographically replicated 
and projected array has the interesting feature 
compared to more traditional optical lattices 
that the position of each trapping site does not 
depend on a relative phase between two inter- 
fering beams. This suggests the potential for 
better long term stability compared to optical 
lattice implementations. 

A particular feature of the Gaussian array 
BoB trap (Fig. [2];) is that a periodic array of 
Gaussian beams creates an array of dark traps 
with the same periodicity, without any extra 



confining walls. This approach would enable 
quasi-magic trapping of 125s atoms on a lat- 
tice with 4 /xm periodicity as detailed in Sec. 
IIVI Conversely the Gaussian interference and 
crossed vortex BoB traps (Fig. [2^,,b) have a 
confining wall around each trap site so that 
there would be two confining walls between each 
site in a replicated array. This implies an ap- 
proximately 50% larger trap period for quasi- 
magic trapping of 125s atoms, which would re- 
duce the number of sites per unit area by more 
than a factor of two. 

As we have shown in Sec.. IIVI quasi- magic 
ground-Rydberg trap shift matching can be 
achieved either by adding a uniform background 
field to the trap designs that have zero inten- 
sity at the trap center (Fig. Ek,b) , or by careful 
choice of the trap parameters of the Gaussian 
lattice trap (Fig. [2J;) which has a finite inten- 
sity at the trap center. These quasi-magic traps 
have no intensity dependent shift for atoms 
at the trap center (or in the motional ground 
state), but do show shifts at finite temperature. 
We have shown in Sec. IIVI that the finite tem- 
perature shifts can be limited to ~ 200 kHz for 
10 /iK Cs atoms, and would be even less for 
colder atoms. Trap shift matching is impor- 
tant for high fidelity Rydberg mediated quan- 
tum gates[5[, furthermore this method may also 
be relevant for high accuracy control of black 
body radiation shifts [22| in optical transition 
atomic clocks. In principle it may be possible to 
improve upon our results by designing a com- 
pensating field with the correct spatial shape 
such that not only the differential shift at trap 
center, but also higher spatial derivatives of the 
differential shift are canceled. We leave this as 
a challenge for future work. 
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Appendix A: Calculation of trap intensity 
distributions 



In this appendix wc document the calcula- 
tion steps used to derive the intensity distribu- 
tions leading to the trapping potentials of Eqs. 
(|l|2|3j) for the Gaussian beam interference BoB, 
crossed vortex BoB, and Gaussian array BoB 
respectively. 



1. Gaussian Interference BoB 

The two Gaussian beams have a phase differ- 
ence of 7T after the Mach-Zehnder interferome- 
ter, and the on-axis intensities are set equal by 
putting P\jw\ = P2/W2 with the power 
and 101 2 the beam waist. The combined inten- 
sity of the BoB trap is 
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2Pi 



wi(z) 



■n(»)] 



W 2 

w 2 (z) 



) e l l kz + k 2R 2 (z) 



where from the properties of TEMoo Gaus- 
sian beams zri,2 = 7ru;f 2 /A, 101,2(2) = 

^-V 1 + ( ^ )2 ' B1 ' 2(Z) = * + % 
^1,2(2) = arct an(^-^) an d p 2 = x 2 + y 2 . Mul- 
tiplying by the polarizability to convert to en- 
ergy units and Taylor expanding about the ori- 



gin gives Eqs. ([T]). 

2. Crossed Vortex BoB 

The intensity of a Lagucrre-Gaussian beam 
can be written as 



where In = 



Ci p — 



and zji = 



2p! 



w(z) 



For the crossed 



vortex BoB we are using, I = 1, p = 0. 

The BoB trap is created by passing two or- 
thogonally polarized beams with separation d 
through a lens of focal length /. After the 
focusing lens, the two beams are rotated by 
±6 = ±arctan(^) in the x, z plane. For not 
too large angles such that we can neglect local 
polarization changes due to the beam focusing 
the BoB intensity is given by 



I(x,y,z) = Ii fi (p + ,z + ) +/i,o(p_,2_) 



with p± = \/ y 1 + (x cos 9 ± z sin 6) 2 and z± = 
z cos =p x sin 8. Multiplying by the polarizabil- 



ity and Taylor expanding about the origin gives 
Eqs. ©. 



Gaussian Lattice 



Each unit cell of the Gaussian beam lattice 
has the same polarization on the upper right 
and lower left corners and an opposite polariza- 
tion on the lower right and upper left corners as 
shown in Fig. &p). We therefore add the fields 
from the beams centered at opposite corners, 
and then add the two intensities. This can be 
written as 



I(x,y,z) 



2P 



\E{x - d/2, y - d/2, z) + E(x + d/2, y + d/2, z)f 



+ \E(x + d/2, y - d/2, z) + E(x - d/2, y + d/2, z)[ 



where each side of the unit cell has length d, normalized field distribution of each beam is 
P is the power of each beam, and the unity 2 2 

E( x , y , z ) = J^e-^e^W-"^. 
w(z) 
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Taylor expanding the potential about the center 
of the unit cell at x = y = gives Eqs. ([3]). 



